Inelastic Confinement-Induced Resonances in Low-Dimensional Quantum Systems 
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A theoretical model is presented describing the confinement-induced resonances observed in the 
recent loss experiment of Haller et al. [Phys. Rev. Lett. 104, 153203 (2010)]. These resonances 
originate from possible molecule formation due to the coupling of center-of-mass and relative motion. 
A corresponding model is verified by ab initio calculations and predicts the resonance positions in 
ID as well as in 2D confinement in agreement with the experiment. This resolves the contradiction 
of the experimental observations to previous theoretical predictions. 



Low-dimensional quantum systems shovif intriguing 
phenomena. For example, a 2D confinement allows for 
the existence of particles with fractional statistics known 
as anyons In ID a gas of impenetrable Bosons, 

the Tonks- Girardeau gas, acquires Fermionic proper- 
ties [2, 3]. Nowadays, one and two-dimensional sys- 
tems can be experimentally realized in trapped ultra- 
cold gases offering a high degree of control 0, 0: Si- 
In low-dimensional systems confinement-induced reso- 
nances (CIRs) have attracted great interest. They are 
universal, since they depend solely on the geometry of 
the trap. In 1998, Olshanii developed a mapping of the 
relative-motion Hamiltonian of a system of two atoms 
confined in a harmonic trap at large anisotropics onto 
the corresponding purely one-dimensional one Q. The 
resulting ID effective interaction strength gio shows a 
divergent behavior at a specific scattering length which 
leads to the formation of a Tonks-Girardeau gas 0, [1] . 
An analogous derivation of the effective 2D interaction 
strength g2D reveals a similar divergent behavior [^. 
The divergences in these models, which are based on the 
relative-motion Hamiltonian in harmonic approximation 
(RMH models), were confirmed in ID by ab initio cal- 
culations of giD and very recently in 2D adopting 
radio-frequency spectroscopy to measure g2D An 
experimental search for the predicted CIRs in terms of 
particle losses and heating was performed by Haller et 
al. Q. In ID, a resonance close to the position predicted 
by the RMH model was found for isotropic transversal 
confinement. However, the experiment considered also 
an anisotropic transversal confinement and found an un- 
expected splitting of the resonance. Furthermore, in 2D 
confinement a resonance was observed at positive values 
of the scattering length. A detailed analysis 12| proved 
formally that these two observations contradict the RMH 
models that predict a single resonance in ID and a reso- 
nance at negative scattering lengths in 2D confinement. 
The immediate question arises what kind of resonances 
were observed in Q. Are the RMH resonances modi- 
fied due to the experimental setup, or are the losses in 
of completely different origin? If so, why are RMH 
resonances not seen in terms of losses? 

In this Letter it will be demonstrated that the reso- 
nances observed in are caused by a coupling of the 



center-of-mass (COM) and relative (REL) motion which 
originates from the anharmonic terms of the trapping op- 
tical lattice. The COM-REL couphng (CRC) and strong 
ID or 2D confinement leads to Feshbach resonances in- 
duced by (avoided) crossings of bound states with COM 
excitation and a state of an unbound atom pair with- 
out COM excitation. Hence, at the avoided crossing a 
molecular bound state can be occupied, since the excess 
binding energy can be transferred to COM excitation 
energy. In view of this energy transfer we refer to the 
COM-REL resonance as an inelastic CIR. The forma- 
tion of molecules is a basic loss mechanism and should 
thus be observable in loss experiments like Q. On the 
contrary, in the spectra of the RMH Hamiltonians [l^ 
no curve crossings appear that could lead to a significant 
occupation of an excited bound state and corresponding 
molecule formation il4| . Because of the absence of an en- 
ergy transfer the CIRs of the RMH models are denoted 
as elastic CIRs in the following. 

The CRC model is introduced, first for ID and later for 
2D confinement, for two identical particles in agreement 
with the experiment. The external potential is assumed 
to be expandable in a power series around the origin and 
separable in the three spatial directions. In COM and 
REL coordinates, r = ri — r2 and R = ^(ri + r2), the 
Hamiltonian is given as 

H{r,R) = TrelW +Tcom(R) + T^relW 

+ VcoM{R) + Uint{r)+W{r,R) (1) 

where Trel and Tcom are the kinetic-energy operators 
of the REL and COM motion, respectively. Vrel and 
VcoM are the separable parts of the potential energy. 
Thus, W contains only the non-separable terms that are 
of the form rfR^ with i g {x,y,z} and n,m G N\{0}. 
^^int (r) is the inter-particle interaction which in the CRC 
model is described by the pseudo potential Uint{r) = 
^^^-^(5(r)^r where a is the 3D s-wave scattering length 
and m the atom mass. In the case of an optical lattice in 
three spatial directions the external potential terms read 



^^REL(r) =2 



(2) 
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(3) 



W^(r,R) = -4 V,^iv?{\kr,)sin^{kR,) (4) 



with = X 



and A the laser wavelength. Vj is the lattice 



depth in direction j S {x,y,z}. In an optical lattice ID 
geometry can be achieved by forming decoupled quasi 
ID tubes which requires a sufficient lattice depth in two 
transversal directions, e.g., x and y. Consequently, it 
suffices to consider the physics of a single tube. 

To estimate the energy spectrum the harmonic approx- 
imation of a single tube with trap frequencies LOj and trap 
lengths dj = yj2h/ (rauij) is considered. The confine- 
ment is characterized by the anisotropics "qx — i^x/^z and 
rjy = ujy/uJz^ which should be sufficiently large. Within 
the harmonic approximation the coupling between COM 
and REL motion vanishes. The wavefunctions are thus 
represented by products of COM and REL eigenstates. 
The harmonic relative motion Hamiltonian possesses a 
single bound state tp^^^r) [l5j. The influence of the 
weak longitudinal confinement on the bound state en- 



ergy E^^^{a) can be neglected [13|. In this case E^^^ 
and a satisfy the relation [l^ 



v/t(l-e-/^*)(l-e-*) 



dt 



(5) 



where e = {Ef^^ - Ea)/{hujy), Eq = f (w^ + ujy), and 
(3 = ujx/^y A general expression for the eigenener- 
gies above the REL motion threshold i?th = f ('^x + 
LOy + ujz) is not known. The states above -Eth are de- 
noted as trap states in the following. The eigenenergy 
^REL ^Yi(^ gj,g|. |.j.g^p g|;g^|;g i^j^^ jjgg i^j^g interval 

[i?th,-Bth + 'i.hujz) such that for the model Ef-^^ is ap- 
proximated by ii^th + [l6| . In the limit as the tem- 
perature r — 0, the higher excitations are frozen out. 
The eigenstates of the COM Hamiltonian factorize as 
$„(R) = (X) (F) (/)„ JZ) with n = {nx,ny,nz) 
and eigenenergies £'C0m _ J2j^x y z + 5)- When 

combining REL and COM motion the energies of the 
bound states become Ef^^{a) + E^'^^ while the energy 
of the lowest trap state is given by Ef^^ + . For 

n 7^ (0,0,0) crossings between the excited bound state 
and the lowest trap state occur for 



^REL 



^REL 



(6) 



where A„ = E^^^ - -^(To'^) is the COM excitation. 
The corresponding scattering length at the crossing is 
obtained from Eq. ([5]). 

In the following, the influence of the anharmonicity of 
the conflning potentials of Eqs. (O, and @ is con- 
sidered. First of all, the non-vanishing term W{r, R) 7^ 



leads to the coupling between REL and COM states of 
equal symmetry. Resonances are only observable if the 
coupling between the crossing states is sufficiently strong. 
Within the CRC model, the matrix elements defining 
the coupling strength between a bound state l^nV" '■'''') 
with COM excitation An and the lowest trap state 
1^(0,0,0) V'l) are W„ = (*n ^'^''WlV'i *(o,o,o))- Since 
VF(r,R) — J2j=x y z^j('''j^ ^j) separates in the spatial 
directions the matrix elements become 



Wn= ('J'nV'^'Wjl*(0,0,0)^l)- 



(7) 



]=x,y,z 



The coupling to a bound state along the weakly confined 
z direction can be neglected, because the bound state has 
an extension db <C dz- Hence, W{z, Z) w W{0, Z) holds 
within the extension of the bound state \z\ < df, and due 
to the orthogonality of the REL eigenstates the longitu- 
dinal matrix element vanishes, (<I'nV'^''^|W^z|'0i'i'(o.o,o)) ~ 
{'^>n\Wz{z = 0,Z)|$(o,o,o))(^^'^IV'i) = 0. With this ap- 
proximation the coupling matrix element ([7]) becomes 



S„^,Q{(bn^4'^''^\Wx\Mo) 



(8) 



As implied by Eq. (jS]), the excitation of the bound state 
A„ must be nonzero for crossings to exist. It follows di- 
rectly from Eq. ([S]) that resonances can only occur for 
excitations of the bound state in a single transversal di- 
rection. For symmetry reasons the excitations must be 
even. The couplings connected to the lowest transversal 
excitations are dominant, because for higher excitations 
the increasing oscillatory behavior of the wavefunctions 
reduces the values of the integrals in Eq. Hence, the 
inelastic CIRs for the quasi ID confinement arise domi- 
nantly from the coupling matrix elements with quantum 
numbers n = (2, 0, 0) and n = (0, 2, 0). 

Based on these selection rules the positions of the in- 
elastic CIR can be predicted. So far, the harmonic ap- 
proximation is used to determine the energy crossings. 
For the excitations of the bound state higher COM states 
are involved so that the influence of the anharmonicity on 
the COM energy cannot be neglected. Therefore, An in 
Eq. ([6]) is corrected within first-order perturbation the- 
ory treating the leading anharmonic (quartic) term of the 
COM optical lattice, — ^ ^j=x ■ 
tion. The resulting energy offset An becomes 



yz ^-Rj, as a perturba- 



j=x,y,z 



(9) 



The CIR positions can now be easily determined us- 
ing Eqs. dS]) and jS]) with n = (2,0,0) and n = (0,2,0). 
Figure [1] shows that the predicted resonance positions 
agree perfectly to the experimental positions of maximal 



3 



atom loss. Note, while in [HI the "edges", i.e. the scat- 
tering length for which the atom loss starts to rise signif- 
icantly, were chosen as resonances positions, in the CRC 
model the conventional positions of maximal particle loss 
are considered The CRC model explains accurately 
the observed splitting of the resonance for anisotropic 
transversal confinement. 
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Figure 1. (color online) Positions of CIRs in terms of the 
scattering length for different values of transversal anisotropy 
in quasi ID. The experimental positions of maximal particle 
loss are compared to the predicted CIR positions of the RMH 
model and the CRC model using the experimental parameters 
(Cs atoms confined in a trap with rjy = 825, Vy = 24.8i7r and 
A = f 064. 49 nm where Er — ft^fc^/2m is the photon recoil 
energy. 

In the following full ab initio calculations of the spec- 
trum of the general Hamiltonian ([1]) are presented to val- 
idate the CRC model. The calculations are performed 
by a full six-dimensional exact-diagonalization approach 
which uses a basis of B splines and spherical harmonics 
[l8| . To incorporate coupling, sextic potentials are used 
which are an accurate representation of a single site of 
an optical lattice The interaction is described by a 
Born-Oppenheimer potential which can be varied to tune 
the scattering length a to arbitrary values [3l ■ While in 
the experiment a quasi ID trap of large anisotropy with 
rjx ~ r]y « 900 is used, the calculations are performed 
with ?7a; ss ?7j, = 10 which is however already well in 
the quasi ID regime (131] . Larger anisotropics would lead 
to a prohibitive computational effort. A fully coupled 
spectrum for transversal isotropy {rix = %) is shown in 
Fig. [2] The complex structure of the energy spectrum 
is highlighting the importance of the derived selection 
rules that specify which of the crossing states couple sig- 
nificantly. The coupling strength is proportional to the 
width of the avoided crossings that can be observed in 
Fig. [3] and the lower graph of Fig. [2l The selection rules 
are clearly verified: significant avoided crossings appear 
only for transversally excited bound states whereas a lon- 
gitudinal excitation leads to negligible coupling. 

While in the case of isotropic confinement, rj^ = rjy, 
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Figure 2. (color online) Spectrum of the full coupled Hamilto- 
nian ([1]) for ^Li atoms confined in a sextic trapping potential 
with Vx = Vy = 35.9Er, rjx = rjy — 10 and A — 1000 nm. In 
the upper part all states bending down to — oo are molecular 
states originating from the REL bound state ■!/)'■''' with differ- 
ent COM excitations. The one in the COM ground state is 
marked red. The magnified part shows the avoided crossing 
responsible for the inelastic CIR which arises from the cross- 
ings of the transversally excited bound states with the ground 
trap state. The kets indicate the dominant contribution to the 
states. Only transversally excited states couple strongly while 
for longitudinal excitation the coupling is very small result- 
ing in almost non-avoided crossings. For isotropic transversal 
confinement only one CIR occurs due to the degeneracy of 
the transversal excitation. 



the transversally excited states are degenerate and only 
a single resonance is visible, for anisotropic confinement, 
Vx 7^ '7j/j this degeneracy is lifted. In accordance with the 
CRC model a splitting of the CIR positions appears in 
Fig. [3] In Fig. |4] the CIR positions of the ab initio calcu- 
lations and the CRC model are shown. The stronger lon- 
gitudinal confinement results in a systematic shift of the 
resonances towards larger values of the scattering length. 
This is mainly caused by the behavior of the energy of 
the first trap state. Although this energy could only be 
estimated within the CRC model the deviations to the ab 
initio results are less than 1%. In contrast, the resonance 
position of the elastic CIR predicted by the RMH model 
is independent of the longitudinal confinement. 

Next, the CRC model is extended to 2D confinement 
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Figure 3. (color online) Spectrum of the full coupled Hamilto- 
nian ([T]) for ^Li atoms confined in a sextic trapping potential 
with Vy = 35. 9S^, rj^ = 11, r/y = 10 and A = 1000 nm. 
The kets indicate the dominant contribution to the states. 
Only transversally excited bound states couple strongly to 
the ground trap state as displayed by the size of the avoided 
crossing. The two CIRs correspond to COM excitations in 
either of the transversal directions. 
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Figure 4. (color online) Positions of inelastic CIRs for dif- 
ferent values of the transversal anisotropy in quasi ID. The 
results of the ab initio calculations are compared to the predic- 
tion of the CRC model using the same parameters ('^Li atoms 
confined in a sextic trapping potential with Vy = 35.9Er, 
rjy — 10 and A — 1000 nm) that, however, differ from the 
ones in Fig. [T] 



which is characterized by a single strong transversal con- 



finement in, e.g., the x direction (i.e. 



While Eq. ([5]) and Eq. © stay valid, the contributions 
of the weakly confined directions to the matrix element 
in Eq. ([7]) can be neglected. Hence, the coupling matrix 
element between an excited bound state and the lowest 
trap state becomes 



(10) 



Only a transversally excited bound state with = 
2, 4, . . . ; TZy = riz = can lead to coupling. Again, only 
the first excitation n = (2,0,0) is dominant, i.e. only a 



single resonance appears. In the experiment [Hj a single 
resonance is observed for 2D confinement indicated by 
maximal particle loss at j- = 0.593 [19]. For the ex- 
perimental trap parameters the CRC model predicts the 
CIR at = 0.595, again, in perfect agreement. This 

ay 

highlights the fundamental difference to the RMH model 
which predicts a resonance for negative values of the scat- 
tering length in 2D confinement. 

As in the ID case the selection rules of the CRC model 
are confirmed by full ab initio calculations. For rj^ = 10 
and rjy — 1 both the ab initio calculations and the CRC 
model predict the CIR position at -j^ — 0.64 showing 

ay 

again that the resonance position depends also on the 
strength of the weak confinement. 

In conclusion, a model describing inelastic CIRs is pre- 
sented which explains qualitatively and quantitatively 
the surprising experimental results of a splitting of the 
resonance under transversal anisotropic ID confinement 
and a resonance at positive scattering length in 2D con- 
finement in the loss experiment [sj . These resonances are 
caused by the coupling to bound states with COM exci- 
tation in the tightly confined direction. At the inelas- 
tic CIR molecules are formed causing enhanced losses. 
The previously known CIRs are elastic and are therefore 
hardly visible in a loss experiment. However, they can be 
observed by measuring the effective interaction strength 
directly ll|. The model has important consequences for 
the experimental realization of a Tonks- Girardeau gas. 
As the inelastic-resonance position depends sensitively on 
the energy of the resonant excited bound state, it is pos- 
sible to prepare a system on the broad elastic CIR where 
(/ID — ^ 00, i.e. the Tonks-Girardeau limit is reached, 
while being off resonant to the excited bound state and 
thus avoiding atom losses Q • In 2D the elastic and inelas- 
tic CIRs are well separated which explains the different 
resonance positions observed in Q and 11 1. Since the 



underlying mechanism of COM-REL coupling leading to 
an inelastic CIRs is of very general nature, we believe 
that they can be observable in various physical systems. 
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